Inertial Effects on Fluid Flow through Disordered Porous Media 
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We study the fluid flow through disordered porous media by numerically solving the complete set 
of the Navier-Stokes equations in a two dimensional lattice with a spatially random distribution of 
solid obstacles (plaquettes). We simulate viscous and non- viscous flow through these idealized pore 
spaces to determine the origin of the deviations from the classical Darcy's law behavior. Due to 
the non-linear contribution of inertia to the transport of momentum at the pore scale, we observe 
a typical departure from Darcy's law at sufficiently high Reynolds numbers. Moreover, we show 
that the classical Forchheimer equation provides a valid phenomenological model to correlate the 
variations of the friction factor of the porous media over a wide range of Reynolds conditions. 
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The application of Darcy's law is the standard ap- 
proach to characterize single phase fluid flow in micro- 
scopically disordered and macroscopically homogeneous 
porous media (U-| . Basically, one simply assumes that a 
global index, the permeability k, relates the average fluid 
velocity V through the pores, with the pressure drop AP 
measured across the system, 
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where L is the length of the sample in the flow direction 
and /i is the viscosity of the fluid. However, in order to 
understand the interplay between porous structure and 
fluid flow, it is necessary to examine local aspects of the 
pore space morphology and relate them with the relevant 
mechanisms of momentum transfer (viscous and inertial 
forces). This has been accomplished in previous stud- 
ies [^|J| where computational simulations based on the 
knowledge of the pore space morphology have been quite 
successful in predicting permeability coefficients of real 
porous materials. 

In spite of its great applicability, the concept of perme- 
ability as a global index for flow, which implies the va- 
lidity of Eq. (1), should be restricted to viscous flow con- 
ditions or, more precisely, to low values of the Reynolds 
number, Re = pVd p / fj,, where p is the density of the fluid 
and d p is the grain diameter. Unlike the sudden transi- 
tion from laminar to turbulent flow in pipes and chan- 
nels, where there is a critical Reynolds number condi- 
tion separating both regimes, experimental studies have 
shown that the passage from the linear (Darcy's law) to 
the nonlinear behavior in flow through porous media is 
more likely to be gradual (see Dullien ji| and references 
therein) . It has then been argued that the contribution of 
inertial forces (convection) to the flow in the pore space 
should also be examined in the framework of the laminar 
flow regime before assuming that fully developed turbu- 
lence effects are already present and relevant to momen- 
tum transport in the system. Recently the transport 



of momentum in two-dimensional porous structures gen- 
erated at the critical percolation point || has been eval- 
uated through the direct solution of the Navier-Stokes 
and continuity equations. It has been shown that, be- 
yond the range of validity of Darcy's law, a nonuniversal 
behaviour should be considered for the critical exponent 
t relating flow permeability coefficients and porosities, 
k oc (e — e c ) , with e c being the critical percolation poros- 
ity. In the present work we show by numerical simulation 
of the Navier-Stokes equations in low porosity structures 
(e > e c ) that the departure from Darcy's law at suf- 
ficiently high Re numbers can be explained in terms of 
the inertial contribution to laminar fluid flow through the 
void space. In other words, we demonstrate that there is 
no need to include turbulence effects to model incipient 
deviations from linearity usually found in permeability 
experiments with various porous materials W . 

Our topological model for the pore connectivity is 
based on the general picture of site percolation disorder. 
Square obstacles are randomly removed from a 64x64 
square lattice until a porous space with a prescribed void 
fraction e is generated. The mathematical description for 
the detailed fluid mechanics in the interstitial pore space 
is based on the assumptions that we have steady state 
flow in isothermal conditions and the fluid is continuum, 
Newtonian and incompressible. Thus, the essence of our 
phenomenological description is the two-dimensional set 
of Navier-Stokes and continuity equations for momentum 
and mass conservation, 
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Here u and v are the components of the velocity vector in 
the x and y directions, respectively. We use the nonslip 
boundary condition at the whole of the solid-fluid inter- 
face. End effects of the flow field established inside the 
pore structure (particularly significant at high Re condi- 
tions) are minimized by attaching a header (inlet) and a 
recovery (outlet) region to two opposite faces arbitrarily 
chosen. At the inlet line, a constant inflow velocity in the 
normal direction to the boundary is specified whereas at 
the exit, the rate of the velocity change is assumed to be 
zero (gradientless boundary condition). Instead of peri- 
odic boundary conditions, we choose to close the remain- 
ing two faces of the system with two additional columns 
of obstacles. This insulating condition reproduces more 
closely the experimental setup usually adopted with real 
rocks and permeameters. For a given realization of the 
pore geometry and a fixed Re, the velocity and pressure 
fields in the fluid phase of the void space and ancillary 
zones are numerically obtained through discretization by 
means of the control volume finite-difference technique 
j?J . Finally, from the area-averaged pressures at the inlet 
and outlet positions, the overall pressure drop across the 
porous realizations can be readily calculated. 




FIG. 1. (a) Contour plot of the stream function ip for low 
Reynolds number conditions {Re = 0.0156). (b) Same as in 
(a), but for a higher Reynolds number {Re = 15.6). 

Figure la shows the contour plot of the stream func- 
tion ip M for a typical realization of a highly porous void 



space (e = 0.9) subjected to low Reynolds conditions, 
Re — 0.0156. In spite of the highly connected pathway 
available for flow at this large porosity value, the predom- 
inant viscous forces in the momentum transport through 
the complex void geometry generates well defined "pref- 
erential channels" of fluid flow. As shown in Fig. lb, the 
situation is quite different at high Reynolds conditions, 
where the degree of channeling is clearly less intense than 
in Fig. la. In the case of Fig. lb, due to the relevant 
contribution of inertial forces (convection) to the flow at 
the pore scale, the distribution of streamlines along the 
direction orthogonal to the main flux becomes more ho- 
mogeneous. This nonlinear effect can be macroscopically 
quantified in terms of the so-called Forchheimer equation 
II, 

- — = a f iV + PpV 2 , (3) 

where the coefficient a should correspond to the recipro- 
cal permeability of the porous material and (3 is usually 
named as an "inertial parameter" . Equation (3) with 
constant a and (3 parameters is not a purely empirical 
expression since it can be derived by an adequate average 
of the Navier-Stokes equation for steady and incompress- 
ible laminar flow of a Newtonian fluid in a rigid porous 
medium Q. Rearranging the Forchheimer equation in 
the form, 

+ (4) 

where / = —AP/Lf3pV 2 and Re' = fipV/ap, we obtain a 
friction factor- Reynolds number type of correlation which 
is presumably "universal" . Indeed, Eq. (4) has been ex- 
tensively and successfully used to correlate experimental 
data from a large variety of porous materials and a broad 
range of flow conditions 0. In Fig. 2, we show the re- 
sults of simulations performed with three realizations of 
the porous structure generated with e = 0.9. After com- 
puting and averaging the overall pressure drops for all 
realizations at different Re numbers, we estimate the co- 
efficients a and (3 and calculate the modified variables 
/ and Re'. In agreement with real flow experiments, 
Eq. (4) also provides a satisfactory fit to all results of 
our computational simulations. Moreover, the point of 
departure from linear (Darcy's law) to nonlinear behavior 
in the range 10 -2 < Re' < 10 _1 of modified Reynolds is 
also consistent with previous experimental observations. 

In an attempt to characterize the influence of inertial 
forces on the flow behavior of a single fluid in highly 
porous structures, we demonstrate here, by direct simu- 
lation of the Navier-Stokes equations, that incipient devi- 
ations from Darcy's law observed in several experiments, 
can be satisfactorily modeled in the laminar regime of 
fluid flow, without including turbulence effects. The re- 
sults of our computational simulations corroborates nu- 
merous experimental data which display a gradual tran- 
sition at high Re from linear to nonlinear flow in the pore 
space. In addition, we show that the physical description 
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underlying the classical Forchhcimer equation provides 
a legitimate correlation for the global friction factor of 
the porous media over a wide range of Reynolds number 
conditions. In summary, whether true turbulence effects 
have been detected or not in real flows through porous 
media, our calculations with the Navier-Stokes equations 
indicate that the Forchhcimer model with constant a and 
(3 parameters remains valid for rather high Re numbers, 
even when convective nonlinearities can significantly af- 
fect the momentum transport at the pore scale, in com- 
parison to viscous forces. These facts might be relevant 
to understand the role played by convection in several 
technical applications involving consolidated and non- 
consolidated porous materials. 
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FIG. 2. Dependence of the generalized friction factor / on 
the modified Reynolds number Re' . The solid line is the best 
fit to the data of the Forchheimer equation. The dashed line 
is the best fit to the data at low Re 1 of Darcy's law. 
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